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FORMULAS FOR THE ELASTIC CONSTANTS OF PLATES WITH INTEGRAL 

WAFFLE-LIKE STIFFENING ‘ 

By Norris F. Dow, Charles Libove, and Ralph E. Hxtbka 


SUMMARY 

Formulas are derived jor the fifteen elastic constants associated 
with bending^ stretching, twisting, and shearing of plates with 
closely spaced integral ribbing in a variety of configurations 
and proportions. In the derivation the plates are considered, 
conceptually, as more uniform orthotropic plates somewhat on 
the order of plywood. The constants, which include the effec- 
tiveness of the ribs for resisting deformations other than bending 
and stretching in their longitudinal directions, are defined in 
terms of four coefficients a, fi, a', and P', and theoretical and 
experimental methods for the evaluation of these coefficients are 
discussed. Four of the more important elastic constants are 
predicted by these formulas and are compared with test results. 
Good correlation is obtained. 

INTRODUCTION 

Growing interest in integrally stiffened construction, 
evidenced by such papers as references 1 and 2 and by the 
large forging press program (ref. 3) and the chemical milling 
process (ref. 4) which will provide facilities for production, 
emphasizes the need for information on the structural char- 
acteristics of integrally stiffened plates. 

A primary requisite for the prediction of structural charac- 
teristics of plates is a knowledge of their elastic constants. 
In the present report, therefore, formulas are derived for 
the fifteen elastic constants associated with the bending, 
stretching, twisting, and shearing of plates with closely 
spaced integral ribs running in one or more directions. The 
ribbing patterns covered by the formulas are illustrated in 
figure 1 and include those considered in reference 5. The 
rib cross section is arbitrary, although special auxiliary 
formulas are given for the rectangular-section rib with 
circular fillets at its base. 

The elastic-constant formulas derived involve four co- 
efficients a, fi, a', p' for each rib which define the effectiveness 
of the rib in resisting deformations other than simple bending 
or stretching in its longitudinal direction. For most purposes 
a reasonably accurate evaluation of these coefficients is 
required. Experimental and theoretical methods of evaluat- 
ing them are discussed. 

As a check on the correctness of the elastic-constant 
formulas, the predictions of the formulas for four of the 
more important elastic constants are compared with experi- 
mental data. 



mj iipu mmp mw 

(d) 

Cb) Longitudinal and transverse, 
(d) Skewed plus longitudinal 
and transverse. 

Figure 1. — Ribbing configurations considered. 

SYMBOLS 

Plane I is defined as the plane in which Nx acts and in 
which €x is measured. Plane II is defined as the plane in 
which Ny acts and in which €y is measured. Plane III is 
defined as the plane in which Nxj, acts and in which jxy is 
measured. These three planes are illustrated in figure 2. 



(c) 

(a) Longitudinal or transverse, 
(c) Skewed. 
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GENERAL SYMBOLS 

coupling elastic constants associated with 

bending and stretching and defined by the 
force-distortion equations (1), (2), (4), and 
(5), lb“^ 

coupling elastic constants associated with 

bending and stretching and defined by the 
force-distortion equations (7), (8), (10), 
and (11), in. 

coupling elastic constant associated with 

twist and shear and defined by the force- 
distortion equations (9) and (12), in. 
bending stiffness in x- and y-directions, 
respectively, in-lb 

twisting stiffnesses relative to x- and y- 
directions, in-lb 

Young'S modulus of material, psi 
extensional stiffnesses in x- and y-direc Lions, 
respectively, Ib/in. 
shear modulus of material, psi 
shear stiffness of plate in xy-plane, Ib/in 
resultant bending-moment intensity in x- and 
y-directions, respectively, lb 


* Supersedes recently declassified NACA RM L53E13a, “Formulas for the Elastic Constants of Plates With Integral Waffle-Like Stiffening" by Norris F. Dow, Charles Libove, and 
Ralph E. Hubka, 1953. 
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■^xv 

N, 

Ny 

N,y 

s 

t 

T 


w 

U,V 

X 

V 

z 


yxy 

€y 

M 

Mx, Ml/ 


M xj M y 
Ml, M2 


hxj hy 

bs 

bs 


bw 

d 


h f h- 

H 

Tw 

liw 

t 

t 

e 


resultant twisting-moment intensity with re- 
gard to x- and y-directions, lb 
intensity of resultant normal force acting in 
2 :-direction in plane I, Ib/in. 
intensity of resultant normal force acting in 
2 /-direction in plane II, Ib/in. 
intensity of resultant shear force acting in x- 
and ^-directions in plane III, Ib/in. 
coordinate, measured parallel to skewed rib, 
in, • 

coordinate, measured perpendicular to skewed 
rib, in. 

coupling elastic constant associated with 
twist and shear and defined by the force- 
distortion equations (3) and (6), lb~^ 
displacement in 2 -direction, in. 
strain energy, in-lb 

coordinate, measured in longitudinal direc- 
tion, in. 

coordinate, measured in transverse direction, 
in, 

coordinate, measured perpendicular to faces 
of skin, in. 

shear strain, with respect to x~ and t/-direc- 
tions, of plane III 

strain of plane I in a;-direction and of plane II 
in y-direction, respectively 
Poisson^s ratio for material 
Poisson^s ratios associated with bending in x- 
and y-directions, respectively, and defined 
by the force-distortion equations (1), (2), 
(7), and (8) 

Poisson^s ratios associated with extension in x- 
and y-directions, respectively, and defined 
by the force-distortion equations (4), (5), 
(10), and (11) 

SYMBOLS REPRESENTING DIMENSIONS 

x-wise and y-wise length, respectively, of 
smallest repeating unit of plate, in. 
spacing of skew ribs, equal to bx/sin 6 or 
by/cos $, in. 

rib spacing (measured between center lines of 
parallel ribs), in. 
rib depth, H—ts, in. 

diameter of largest circle that can be inscribed 
in cross section at intersection of rib and 
skin, in. 

distance from planes of zero strain to rib 
centroids, in. 

overall height of rib plus skin, in. 
radius of fillet, in. 
corner radius, in. 
thickness, in. 

average or equivalent thickness, in. 
angle of skewed ribbing, measured from the 
longitudinal direction, deg 


SYMBOLS USED IN EQUATIONS FOR ELASTIC CONSTANTS 


a 

constant used in equations for calcu- 
lating a'uL 


cross-sectional area (including fillets) 
of aj-wise, y-wise, and skewed ribs 
(Aw^ includes area of two ribs), 
sq in. 

Aw 

general symbol for Aw.^, Aw^, or Aw^ 

y, g, ^ 

constants used in equations for calcu- 
lating auL and 


cross-sectional moment of inertia of 
x-wise, y-wise, or skewed ribs 
about their centroids {Iw^ is twice 
the moment of inertia of a single 
skew rib), in.'^ 

ki, kiiy kiu 

dimensionless distance from middle 
surface of sheet to planes I, II, 
and III, respectively, expressed as 


fractions of the overall height H 

kw^y kw^) kw^ 

dimensionless distance from middle 
surface of sheet to centroid of x- 
wise, y-wise, or skewed rib, ex- 
pressed as a fraction of the overall 
height H 

{ constants used to locate the effective 

OtLL) <^exp} OCxjL 

centroid of a rib for resisting bend- 

Otx) ^y) 

1 ing in its transverse direction 

a 

general symbol representing a^y a^, 
or a, 

ot\, a! y, ol\ 

constants used to locate the effective 
centroid of a rib for resisting 
twisting 

a' 

&LL, ^exp, ffvL 

general symbol representing a'x, 
or a\ 

[ constants used to define effectiveness 

of a rib in resisting stretching in 

its transverse direction 

& 

general symbol representing &yi 

or 

fi'v, H'a 

constants used to define effectiveness 
of a rib in resisting shearing 

/S' 

general symbol representing P'vi 

or 

SUBSCRIPTS 

L 

longitudinal 

S 

sheet or skin 

T 

transverse 

W 

rib (web) 

s, X, y 

indicate application to skewed, x- 
wise, or y-wise ribs or directions 

LL 

lower limit 

UL 

upper limit 

exp 

experimental 
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DEFINITION OF ELASTIC CONSTANTS 


If the rib spacings of integrally stiffened plates such as 
those shown in figure 1 are small in comparison with the 
plate width and length, it is plausible, for purposes of study- 
ing overall or average behavior, to assume that the actual 
plate may be replaced by an equivalent uniform orthotropic 
plate. Figure 2 shows an infinitesimal element of the 
equivalent plate subjected to bending moments of intensity 
Mx and My, twisting moments of intensity M^y, stretching 
forces of intensity Nx and Ny acting in planes I and II, 
respectively, and shearing forces of intensity Nxy in plane 
III. The locations of planes I, II, and III are arbitrary. 

The behavior of the element can be described by a set of 
force-distortion relationships in which elastic constants 
appear. Such relationships for special rectangular ortho- 
tropic plates having their axes of principal stiffness parallel 
and perpendicular to their edges, as considered herein, are 
obtainable from reference 6. If deflections due to depth- 
wise shear are assumed to be negligible as is customary in 
ordinary plate theory, the following equations (eqs. (1') 
to (6') of ref. 6) are obtained: 


Mxy I rp■f^T 

dxi>y ^ 


€y=-C,yM,-CyyMy-f^^ 


y^=2TMry+ 


Gxv 


( 1 ) 

( 2 ) 

(3) 

(4) 
(5j 
( 6 ) 


where and are the curvatures, is the twist, 

€x and €y are the extensional strains in planes I and II, re- 
spectively, and is the shear strain in plane III. 

According to these equations, fifteen constants are needed 
to establish the force-distortion relationships — namely, two 
bending stiffnesses Dx and Dy, a twisting stiffness Dxy, two 
stretching moduli Ex and Ey, a shearing modulus Gxy, two 
Poisson's ratios Mz and fXy associated with bending, two 
Poisson's ratios •/x and ix'y associated with stretching, four 


CUUpiiiig tei Ills ^ ^ xyy Eyxj 


id Cyy associated with be^din: 


and stretching, and one coupling term T associated with 
twisting and shear. Not all these constants are independent, 
however; for example, as a consequence of the reciprocity 
theorem for elastic structures, fiy=DyfjLxlDx and tx\=Ey^i'xlEx, 
The form in which the force-distortion relationships have 
just been given is not the most convenient form for some 
applications, particularly for buckling calculations. For 
such purposes a more suitable form is obtamcd when the 


z,w 



plane I 

Plone n 

Pione in 

Figure 2. — Forces and moments acting on element. 

first three equations are solved simultaneously for Mx, My, 
and Mxy and these expressions are then used to eliminate 
Mx, Myy and Mxy in the last three equations. The six new 
force-distortion equations thus obtained are 


My=-D, (^+i^y^^+CnN,+C,,Ny 

«x-On + 


bx? 




yxy — 20* 


Er 

Nry 


Ei 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


where Hy^DiHxjDi and ix 2 =EiHilEi. 

Of the fifteen elastic constants appearing in equations (7) 
to (12), two, Hx and Hy, were also in the original set of force- 
distoi tiou equations. The remaining constants (O., D., D*, 
E^, Ei, Ok, Ml, Mi, 0„, Oi 2 , O 21 , O 22 , and O*) are “new.” The 
algebraic relationships between the new and the original 
elastic constants are given in appendix A. 

METHOD OF ANALYSIS 

The anal 3 'sis is made for a plate with the general pattern of 
ribbing shown in figure 3 (a), which includes, as special cases, 
the patterns of figure 1. A typical repeating element of the 
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plate is indicated by the short-dashed rectangle in figure 
3 (a) and is shown three-dim ensionaUy in figure 3 (b). 

The analysis is based on the assumption that each of the 
four rib segments shown in figure 3 (b) may be replaced by 
three orthotropic sheets of material parallel to the skin, each 
one covering the entire area b^by and each fastened to the 
skin by means of many hypothetical, perfectly rigid, infinites- 
imally small bars imbedded perpendicularly through the skin 
and sheets (see fig. 4). The substitute sheets are assumed to 
offer no interference to one another. (The rib is understood 
to include any fillet material but no part of the skin.) The 
properties of the three substitute sheets are so chosen that 



Figure 3. — Repeating element of plate with integral, waffle-like 

stiffening. 


one sheet (labeled © in fig. 4) represents only the effective- 
ness of the rib in resisting stretching and bending in its 
longitudinal direction, another (labeled ®) represents only 
the effectiveness of the rib in resisting stretching and bending 
in its transverse direction, and the third (labeled ©) repre- 
sents only the effectiveness of the rib in resisting shearing and 
twisting relative to its longitudinal and transverse directions. 
(The transverse direction, as used herein, is the direction in 
which tw is measured, see fig. 3.) In order for the three 
substitute sheets to accomplish their purpose, they are 
assigned the following properties: 

(a) Sheet © has a volume equal to that of the rib segment 
it replaces, with its center of gravity at the same level as 
that of the rib. Its stretching or compressing modulus of 
elasticity in the direction of the rib is E and its modulus 
transverse to the rib is zero. Its stiffness per unit width for 
bending in the direction of the rib is equal to the bending 
stiffness of the rib about its centroid divided by the rib 
spacing (i.e., bx for a y-wise rib, by for an x-wise rib, and 6, 
for a skew rib, fig. 3(a)), whereas its bending stiffness in the 
direction transverse to the rib is zero. The shearing and 
twisting stiffnesses and Poisson^s ratios of the sheet are 
assumed to be zero. 

(b) Sheet ® has a volume equal to some fraction of 
the volume of the rib segment, with its center of gravity at 
some distance aH above the middle surface of the skin. 
The modulus of elasticity for stretching or compressing in 
the direction transverse to the rib is whereas that in the 
longitudinal direction of the rib is zero. The bending, shear- 
ing, and twisting stiffnesses, and Poisson’s ratios for sheet © 
are all assumed to be zero. 

(c) Sheet © has a volume equal to some fraction of 
the volume of the rib segment, with its center of gravity at 
some distance a'H above the middle surface of the skin. 
Its modulus of elasticity for shearing relative to the longi- 
tudinal and transverse directions of the rib is whereas its 
twisting stiffness relative to these two directions is zero, 
as are the stretching and bending stiffnesses and Poisson’s 
ratios. 

On the basis of the foregoing assumptions, the integrally 
stiffened plate has been converted to a more homogeneous 
plate somewhat on the order of plywood. The assumption 
of rigid bars connecting the substitute sheets and the skin 
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is equivalent to the assumption that material lines normal 
to the surface of the plate before deformation remain straight 
during deformation. If it is further assumed that these lines 
remain perpendicular to the surface of the plate and that the 
stresses are in the elastic range, any of the methods used for 
ordinary isotropic plate analysis may be readily extended to 
the present idealized structure. 

For the present purpose an energy method is adopted to 

determine the six forces and moments necessary to nmintain 

, ., j .. , . .. b^w b^w 

the prescribed uniform deformations 

and The equations obtained for these forces and mo- 
ments in terms of the distortions are put in the form of equa- 
tions (1) to (6) to yield formulas for the original elastic 
constants or in the form of equations (7) to (12) to yield 
formulas for the new elastic constants. 

The details of the analysis and the derivation of the elastic 
constants are presented in appendix B. The formulas 
obtained for these constants are presented in the following 
section and the evaluation of a, ct , and is discussed in 
two succeeding sections. 


FORMULAS FOR ELASTIC CONSTANTS 

In this section the formulas are presented for the calcu- 
lation of the fifteen elastic constants appearing in equations 
(1) to (6) and the thirteen new constants appearing in equa- 
tions (7) to (12). The formulas are presented for the most 
general type of plate considered, which is illustrated in 
figure 3. For plates with one or more sets of ribs omitted, 
the formulas also apply when the terms representing the 
areas and moments of inertia of the omitted ribs are set 
equal to zero. 

The formulas for the constants in the original force- 
distortion equations (1) to (6) are as follows: 




D.=Em 


E,=EH 


6i 


,=eh\ 


Axylx\ 


\jxu “h 4-4x1/ 


-Arm)*] 


(13) 


Dy = EEP^Iy-^ (14) 

D^=EW(^ (15) 


/ A,\IJy-I,^)-AMxAy{h-hyQCy-hy-A,^\_^xhih-hf^-2^^^I,(JCy-h)QCv-h)^-AyL{ky-hy\ 
Ay{LIy-U)+AyAyIy{h-hy-A.{k,-kl) [A,Iy(Jc-kl) -2AyI.{TCy-k,) +A,Ay{k~kd (^,-^.)*] J 

I A.^ {IJy-I,^ -AMyAy(k-k,y(JCy-k.y-A.^ ^^Iy(jCy-k.y +2^^ I,{h-k,) {ky~k,) +Ayly(ky-k,y~\ \ 

1 Ay(IJy-I,^+A,Ayh(fiy-ki,y-A.Oi.-ht) [A,Uk,-k„) -2AJ,ih-k.) +A,A,(k,-kii) (h-k.y] J 


(18) 


Mx= 




IyA,^—A,^Ay(jCy—kyy 


(19) 






I,A^—A^Ay(JCy—k^* 


J A,(I,Iy-I,^+A^^Iy(h-ki) (h-k,)+A.AyI^(ky-k„)(ky-k.)+AyAyI.(ky-ki)(ky-k„)- \ 
^ I -4,'^/,(X, — ki)(jc, ki^ -^A,AxAyQcz k-^(jcy kji)(jcx k^(ky k,) J 

" Ay(IxIy-I.^+AxAyIy0ix-kiy-A,(Jc,-ki)\A,Iy{k.-h)-2AyI,(Jcy-k.)+A.Ay{k,-ki)(jcy-k.y] 


( 20 ) 


( 21 ) 


f A.{IxIy-I.^)+A.AxIy (kx-h) (kx-ic.) +A,AyIx(jcy-ku) (h-k.) +AxAyI.(kx-h) (ky-kn) - 


j * — y — • / I a \ > * ir •“ N If ■ 

[ A,u,(k.-h) (i,-An) +A.AxAy(kx-kj) (^y-k„) (jiy-k.) 
I* » T . . _ /T ^ \2 . /V i \ r - r #T ~i \ 7 T~a r /T 


'Ax(IxIy-I.^+AxAyIx(ky-kny-A.(,h-ku)[A.Ixik,-kn)- 2 AxI.ikx-h)+A.Ax(k,-kn)(kx-Q ] 


( 22 ) 


o*= 


EIP 


, AxAykx—A,% A,Ay(ky—k,) 

31 J2 


(23) 
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EH^ 






, AxAyky As^k^ 

^II -FT 




(24) 




1 

EH^ 


A,Ay{ky —k, 

A, 


',) Z',. A^Ayh-A,'‘k\ 

"“"r — 27 — ) 




A, 


(25) 


r 

^VV-EH^ 


j AjAyky A^ks A^Axikx k^ 
Tl f^y T o 


A, 


As^Ay U { 


rp 1 P ^IIl) 

TTfT-TO 




] 


(26) 


(27) 


The formulas for the constants in the new equations (eqs. 
(7) to (12)) are as follows; 



(28) 


(29) 

D,=Ew(^-f^ 

(30) 

i+h 

X 

(*q 

II 

(31) 

/a 2\ 
E2^EH(^) 

\Ax / 

(32) 

G,=EH{Axy) 

(33) 

A, 

Ay 

(34) 

A, 

(35) 

Cn^H^ki 

(36) 

1 1 

1 

IXI 

II 

(37) 


(38) 

C22 -H [kn 

(39) 


Ck—H {kxy A^iii) 


( 40 ) 


where 

E Young's modulus of material, psi 
H overall height of skin plus ribs, in. 

The quantities As^ /,, Axj Ay^ and Axy, kxy ky, and 
kxyy Ixy lyy I sy aud I xy appcaring in equations (13) to (40) are 
defined by the following equations: 


ij, — 


(41) 

// = 7,^/ +A,AxAyikx~k, )(ky~ks) (42 ) 


cos“ 

1 is , ^ Awjh. 


(43) 


sin^d cos^ d 

1 + M 


V Jhx Awjhg / 

^+— ^-(sin^0+/3. cos^e+ 


) 


(44) 


^sin^0 cos^d+fis sin^0 cos^0— 


3' 

M+m' 


») 


(45) 


A 


1 is I o/ 1 Awjby 1 Afyjbx 

zj\P a: 0/1 I A rr ■ r P I 


^ 2(1+m) i? ' ^"2(1 +m) H ^'"■'2(1+m) H 

^sin^e cQs^6+fi, sin^e cos^e+/3'j cos^20j 


(46) 


h=l (k.x) ^ 


2 

^socs sin^^+/S',a'« —a — sin^0 cos^0 

1+M 


')] 


(47) 


7 1 ffl , ^ I ^Wylbz ,J V Awjb, /. 


2 

/SsCKa cos^^+/3's«', ~ — sin^0 cos^0 
1+M 


')] 


fcwr, sin^0+ 
(48) 


— 1 ^ Ay^ Jbi / — 

4=^ L cos^e+jSjaj sin^'e cos^S— 


’ set's - r- sin^e cos^d 
1+yu 


)] 


(49) 


1 r 1 A\v Jhy 1 /^x 

_V /Q' - ! 


^ »■ o7l _L .A ‘t/ (“ V O/'l 'l A (^4 + 


2(1 +m) H 


'2(1+/.) 


sin^e cos’^fl+lSsas sin^e cos^^+ 


H 


2 ( 1 +^) 


]} 


(50) 
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r_ 1 1 , 1 ts rr\i 

12(1-^") (a) + tp H> “>* *+!_„• a W + 


(I, -fe)*+|S,^- 

r/r r 




fo)* cos^e-|-ft(o,— Ji)* sin‘ 0-f 


ff'tia't—kxY 8in*d CO8*0^J 


12(1 


1 is /T 


( 1 ,- 4 ,)’+ 

4»)* sin*(»+/S,(a,— 4,)* cosV+ 
fi'$(a',—\y COS^e^J 


( 51 ) 

(^k)*+ 


(52) 


/.= 


M /^sV 

i-M*) Vh; 


V . 


^ /'I 




12(1-m*) \HJ ^ m H 

co8‘*e+i8,(a,— t,)' 8in*e C08^e— 
|8'. {a',—k,y (^Y^ sin*^ cos*e^J (53) 

6(1+m) \HJ 


tsV , . W** 


+4 




sin*^ C08*S- 


1+1^ H 


(kj^+ 


Q> 2 AwJbj, , y^S' 2 ■A.yrjhx , ■vj 

^ * l-f-/i JE? V 1~1“M H ^ ^ 

AwJb.C- 


H 


■{^ (kw,—k„y si 


sin^e cos^e+ 


i3,(a,— cos^^+/3^(a^~0^|^ 2(l^ '^ cos^20j 

(54) 

where 

bxj by, bg the spacing of the x-wise, y-wise, and skew ribs, 
respectively, in. 

0 the angle of skew of the ribbing, deg 

H the overall height of skin plus ribs, in. 

ts the thickness of the skin, in. 

fi Poisson^s ratio for material 

Equations •(43) to (54) contain the quantities A^^, 

and Awgf ^Wyy ^Wgt and Iwyy and Iwg which define 
the areas, iocaiiuns of ceatroids, and moment ? of inertia of 
the ribs. For rectangular ribs with circular fillets, as shown 
in figure 4, these quantities are given by the equations 




(Eq. (57) contains a factor 2 to account for the fact that 
there are two ribs in the skewed direction — one at an angle 
+0 to the x-direction and the other at an angle — ^ to the 
x-direction.) 

H (58) 

(59) 






H 


is 

H 

(60) 




Iwjhi 

~w~ 




+ 




(63) 


The values of ki, ku, and km depend upon the locations of 
the centroids of the forces Nx, Ny, and Nxy, respectively, im- 
posed upon the plate element. (See fig. 2.) For the impor- 
tant case in which Nx acts in such a plane that it produces 


no curvature and Ny acts in such a plane that it produces 

W 

2) and, therefore, 


no curvature Cxx and Cyy must equal zero (see eqs. 1 and 


f_ AxAykx A^kg-A ^^xAgAyiJcy k^) 

7? 


(64) 


533978 0 -60 -2 
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1 Ax-^yky .^s S-^x(J^X 

— — 

A 2 


(65) 


Similarly, for the case in which Nxy acts in such a plane that 
it produces no twist T must equal zero and, therefore, 


kui — iCxy 


( 66 ) 


If Nx and Ny do act in such planes that they produce 

curvatures ^ and the actual locations of the forces 

(planes I and II) must be known if constants (such as Ex^ 
At'z, etc.) which depend upon the locations of the applied 
forces are to be evaluated. 


The average or equivalent thickness of integrally stiffened 
plates having rectangular ribs with circular fillets in the 
various configurations considered herein may be calculated 
from the following formulas: For simple longitudinal or 
transverse ribbing (fig. 1(a)), 


H 


= 1 - 






(67) 


or 








(68) 


For combined longitudinal and transverse ribbing (fig. 1(b)) having a corner radius Rw^^w and with blending 
smoothly with at the corners, 

rf ^ TT 1 7 


ts 


(69) 


For the special square pattern of longitudinal and transverse ribbing having tw^=tw and rw^=Vw , equation (69) 
reduces to 


H H /b.V 


¥.(h 


) 


(69a) 


For skewed ribbing (fig. 1(c)), again for Rw^'T'w, 


H 




(70) 


For 45° skewed ribbing, equation (70) reduces to the form equivalent to (69a); thus, 

1= , 


H' 


-(^Y 

\^s) 


(70a) 


For combined skewed ribbing and transverse (or longitudinal) ribbing in a pattern of triangles, again for Rw'>Twj 




t (uXH-u ») - ») - *' 

T(ir 


+ 


CSC 26 






( 71 ) 
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For the special equilateral triangle pattern having tw=tw^=tw and tw=Tw=Tw, as considered in reference 7, equation 
(71) reduces to 


I . (I -‘Id-'-' d)’{® [d“‘= ?]} 

R-^ R /6,V 


^5 \^S/ 


(71a) 


For combined longitudinal and transverse and skewed ribbing as illustrated in figure 3(a), for Rw'>Tw as before, 

tan«)]+2{2,-[(l+,ec *)>+(l+csc | 




is ts ts 


-0.429 1 ( 1 CSC cot «)+( '-p )\ «c tan ») 

^ (?■') (I I" §+°°‘{^}])} 


Hhh 

fs 


cot|+(^) cot(5^)-h(^*)(cot<?+eot{90-<?}+cot|-^cot{^})- 


(72) 


HhK 

and, finally, for the special case of combined longitudinal and transverse and skewed ribbing illustrated in figure 1(d) 
having tw^—tw^=tw^f and 6x=6y=L4146|, equation (72) reduces to 


R~ H /kV 


\^sJ 


(72a) 


EVALUATION OF a AND /3^ 

EXPERIMENTAL EVALUATION 

The coefBcients a, jS, a', and /S' occurring in the equations 
for the elastic constants express the effectiveness of a rib for 
resisting deformations other than bending and stretching in 
its longitudinal direction. For the evaluation of a and /3 for 
a given set of ribs (longitudinal, transverse, or skewed) prob- 
ably sufficient accuracy will be achieved from a direct experi- 
mental measurement with a simple model having one set of 
ribs with cross section and spacing that duplicate those of 
the ribs for which the coefficients a and ^ are being sought 
and with a value of ts equal to that of the actual plate. 

A double specimen of the type shown on the right-hand 
side of figure 5 may first be used to evaluate through a 
tension test and, then, one-half of the specimen may be used 
to evaluate a through a bending test, as illustrated on the 
left-hand side of figure 5. The use of a double specimen for 
the stretching test is suggested because the symmetry will 
eliminate localized bending of the skin between ribs and facil- 


tSubsequent work has shown that the approximation is adequate for all practical 

purposes (see ref. 8) . 


itate the measurement of overall strain. Because of the pre- 
vention of localized bending, the value of ^ should be some- 
what higher than that which would be obtained by stretching 
a single specimen like the one on the left-hand side of figure 5. 
However, such an overestimate of may be desirable if the 
actual plate has ribs in more than one direction, because then 
the localized curvatures associated with one set of ribs will 
tend to be reduced by the presence of the other ribs. 

The length-width ratio of the specimen should be great 
enough so that any end grips or heavy end sections will offer 
negligible resistance to transverse contraction in the stretch- 
ing lest and to the development of transverse curvature in 
the bending test. Furthermore, the width of the specimen 
should be sufficiently large compared with the rib spacing so 
that the percentage of the specimen subject to shear-lag 
effects arising at the rib ends is small. 

The use of these tests for the evaluation of a and ^ are now 
described in detail. For ease in discussion, the ribs for which 
a and jS are being sought are assumed to be oriented in the 
^/-direction as shown in figure 5. After the values of and 
Py have been determined, however, the subscript y should 
be changed to x or s if, in the actual plate, the ribs under 
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right-hand side of figure 5 are -^— ^=Ny=0. Substi- 
tuting these conditions in equation (10) and making use of 
equations (31), (41), (43), (44), and (45) gives 


A * 

=EH^ 

Ay 

Ay 


=EH 


(i^ 


1 1 


H 


0x\f 1 

Al-M^g'*' H 

1 h I 

H 




Solving for /3» gives 


Ny tf, 




is I 

H 


where, for rectangular ribs with circular fillets, — ^ — is as 

given by equation (56). If the value obtained in the stretch- 
ing test is used for NxjEHtx in the right-hand side of equation 
(74), an experimental value of or is obtained {tx 

is the x-wise strain averaged over at least one multiple of 6x)- 
The conditions of the bending test illustrated in the left- 
hand side of figure 5 are Nx=Ny=My—0. Substituting these 
conditions inequation (1) and making use of equations (13), 
(19), (42), (47), (49), (51), and (53) gives 


Mx_ Mx l__ 

5?- d;- mp _^ 

* ^ 2 


k.y ^*(^ 2) 


where 


I — ^ I ^ {Ic V I fl )2 

h=^R 

Ay H 


IyA,^-A,Ayky 


I,^=I,A,^AA,AyA^JCy 
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Solving for gives 



where, as before, for rectangular ribs with circular fillets, 
— ^ — is as given by equation (56). 


Substituting for the value obtained in the bend- 

ing test, and for /3, the value obtained from equation (74) 
permits equation (77) to yield an experimental value of 

cty is the a:-wise curvature averaged over at least one 

multiple of 6,^. The quantities Z,*, A,, Ay, A„ Icy, /„ /, 

are obtained from equations (41), (43), (44), (45), (48), (52), 
and (53), respectively, with Aw^=Aw=Iw,=Iw,—0-, thus. 


mens shown in figure 5 for their small deformations under the 
assumption that plane sections perpendicular to the skin and 
perpendicular or parallel to the direction of ribbing remain 
plane. The results of such an analysis of the two situations 
illustrated in figure 5 are as follows: 

For the double specimen on the right-hand side of figure 5, 


Ny 

EHu 


( 1 — 


tsIH 


(80) 


For the single specimen on the left-hand side of figure 5, 


A,*=AyAy-A,'‘ 


A _A. Ajry/bx _ 

T TJ 

fi £1 






- _ 1 Awjhy - 
Ay H 


(78) 


I A ly^Jhy ^ ^ I 

^ ^*“12(1 -A.*) vs; J 

where kwy 'i^ as given by equation (59). 

THEORETICAL EVALUATION 


My 1 

«'-<‘V(^)(f)*+^ »1) 

where /, g, and j are geometric properties of segments of 
length hx of the cross sections shown. The symbol I rep- 
resents the moment of inertia of such a segment about its 
centroid, g is the integral, taken in the a;-direction, of the 
reciprocal of the local thickness measured in the 2 -direction, 
and/ is ts^ times a similar integral of the cube of the reciprocal 
of the local thickness. When the ribs are rectangular with 
circular fillets, these quantities are given by the following 
formulas: 


Accurate theoretical analysis of the situations depicted in 
figure 5 is difficult. However, it is possible to obtain values 
of a and that underestimate or overestimate the stiffness 
of the specimens. An underestimate is obviously obtained 
by assuming no part of the rib to be effective in resisting 
transverse stretching or bending in a direction transverse to 
itself. A lower limit-value of ^ is, therefore, 

(79) 

Wlicu ^ is taken as zero, the value of a is immaterial. 

An overestimate is obtained by analyzing the two speci- 


FP^12\hJ^H 


1 ^ ^ (^wjh 

o tr“r 4 \ 


2H^ts\ H 




Iwjbx AwJU 


TP 


H 


2^+^" 


H(AwJb, 

ts\ H 


) 


+ 


+ 1 




( 82 ) 
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^ ^ ts^ ts 




(83) 


-h.JZi-0 

~ts ts ^ ts^ts\Hj^J 


f- 


where g' and f are functions of the ratio of fillet radius to 
skin thickness given by the equations 


(/'=4 


r- 


(85) 




f 

— “ 


h'~ 2 . 

1+2 

/ 1 

(84) 


V-(r+2) 







1 


-^+4^+4 

Vw 


2+ 


rw\ 


Tw 




tan“* -y/ 


1+2 




( 86 ) 


and 


Iwjhx 


is as given by equation (62). 


b^w 


The values of NxlEH^x and MJEIP obtained from 

equations (80) and (81) may be thought of as experimental 
results and they may therefore be substituted in equations 
(74) and (77) to obtain values of and auz, corresponding to 
an overestimate of the stiffness of the specimen. 

A lower overestimate of stiffness can be obtained by analyz- 
ing, on the basis that plane sections remain plane, the single 
specimen on the left-hand side of figure 5 for both N^IEHex 

and MxjEIP ^ and thus including the localized bending 

that occurs during stretching. Besides being more conserva- 
tive, the resulting values of auL and Pul would also be more 
appropriate if, in the actual plate imder consideration, there 
were really only one set of ribs. An upper-limit analyBis 
conducted entirely on the specimen on the left-hand side of 
figure 5 would yield the following expression to be used in 
place of equation (80) : 


Nx 


ts/H 


EHe 


12 (1 


2,(9 


(87) 


Awjhz 

H 


(!) 


where h is ts times the integral, taken over a length hx in the 
x-direction, of the square of the reciprocal of the local thick- 
ness. For circular-filleted rectangular-section ribbing, 


^ t, u f. t, 




( 88 ) 


where W is given by the equation : 


tan" 


‘ V‘+^v 


y (89) 


b^w 


Equation (81) would stUl be used for MJEIP 
EVALUATION OF a' AND /3 t 

The coefficients a' and which define the effectiveness of 
a rib in resisting twisting and shearing relative to its longi- 
tudinal and transverse directions, are not as readily measured 
experimentally nor as readily bounded by an upper limit as 
a and P, although, of course, a lower-limit stiffness is ob- 
tained by equating /S' to zero. 

An approximate evaluation of a' and /S' may be made by 
assuming that the same volume of rib material resists shear 
as resists transverse stretching, that is, 


P\ = Pv 


(90) 


and then by determining from computations where this 
material must be placed (as measured by a') in order to give 
the proper torsional stiffness as determined with the aid of 
reference 9. The computation of a' is now described in 
detail. 

Consider an element, like the one on the left-hand side of 
figure 5, having only y^wise ribbing and subjected to a pure 
Mjj, loading. From equations (3), (15), and (54) the value of 
a\ can be obtained in terms of the measured or computed 

ratio Mxyl ^ as follows: 


bx by 


Mxy _ 

b^v) 
bx by 


-Dr 


-^EHUxy 




P y 1 .. TJ V 


l+M H 


a-' y — T(ty)‘‘ 


where 


Awjhz 

P y jJ y) 


(91) 


(92) 


Solving for aj gives 


H\r 


2\^^Aw^jhx ts 


H 


H 


Mz, 




■m 


2(1 +m) bx by 


(93) 


tA comprehensive evaluation of a' and /8' is now available in reference 8. 
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The value of the ratio M- 




to be inserted in 


2(H-/j) dz by 

equation (93) can, in the absence of test data, be derived by 
an adaptation of the method used in reference 9 for com- 
puting the torsional stiffness of I-beams and H-beams, which 
gives 




EH^ b^w 


2 (1 -1 -m) bx by 


(94) 

where d is the diameter of the largest circle which can be 
inscribed in the cross section at the junction of the rib and 
skin and can be computed from the formula 


d_ 

ts 




(95) 


The constant a in the last term of equation (94) depends on 
twjts and rwjts- The value of a is obtainable from figure 7 

{ '^w \ 

of reference 9 or, when -7^ >0.61—0.23 from the 

Is \ ts / 

following formula; 


a=0.094-(-0.070 


rw, 

ts 


(96) 


The meanings of the various terms within the parentheses 
of equation (94) are apparent: represents the con- 

tribution of the skin, considered as an infinite plate, to the 

twisting stiffness of the waflBe; | (M) Cf) 

is similarly representative of the twisting stiffness of the rib; 

/tw V 

the term with — 0.105( corrects for the fact that the rib 
is actually not infinitely deep; and the term with o( 


represents the additional stiffness due to the fillets. The 

value 0.105 is based on the assumption that - — ^>2.3; for 

values of 2bwyjtw^ less than 2.3, the number '0.105 should be 
replaced by the number obtainable in figure 3 of reference 9 
with the abscissa label b/n replaced by the label 2bwyltwy^ 


MEASURED VALUES OF ELASTIC CONSTANTS 


As a partial checkon the theory, experimental measurements 
were made of the stretching stiffness £*1, bending stiffness 
Dxy shearing stiffness Gky and twisting stiffness D^y of plates 
with integral ribs running either longitudinally or trans- 


versely (fig. 1(a)) or skewed (fig. 1 (c)). The procedures 
used for the measurement of and were essentially the 
same as those described in reference 6 for sandwich plates. 
The measurements of Ei and 6rjt were made with long-gage- 
length resistance-type wire strain gages mounted in the four 
corners (or diagonally on the four sides) of square-tube com- 
pression or torsion specimens similar to the square tubes of 
reference 5. The compression specimens were tested in the 
1,200,000-pound-capacity testing machine and the torsion 
specimens in the combined load testing machine of the 
Langley structures research laboratory. 

The experimental values obtained for the stiffnesses are 
indicated by the circles in figures 6 and 7. In figure 6 the 
stiffnesses are plotted against the angle of skew of the ribbing 
(with 0=0® and 0=90® corresponding to purely longitudinal 
and purely transverse ribbing, respectively) for pla jes having 
nominally the same weight. In figure 7, for a given angle of 
skew (0=45®), the variation of the elastic constants with 
skin thickness is plotted. The relatively large scatter in the 
test data is due to the fact that the plates used were sand 
castings and, hence, had appreciable variations in thicknesses 
from one specimen to another and also within each specimen. 

For comparison, theoretical values of the four elastic 
constants were computed from equations (31), (13), (33), and 



Figure 6. — Calculated and experimentally measured elastic constants 
for plates with integral, waffle-like stiffening skewed at angles of 
± 0 to the longitudinal direction, having 6 jf= 0.2 in., E= 10.7 X 10® ksi, 

ju=0. 32, and having the following proportions: ~=4, ^=2, — =2; 

ts tw 

in addition, for ^=0® or 90®, ^=0.4 and for O®<0<9O®, ^=0 2 

Os Os 
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Figure 7. — Calculated and experimentally measured elastic constants 

for plates having integral, waffle-like stiffening skewed at angles of 

±45° to the longitudinal direction and having ^=2, ——2, ^=0,2, 

tw Tw Os 

6h^= 0.2 in., .E=10.7X103 ksi, and m = 0.32. 

(15) and are plotted in figures 6 and 7. The lowest curve in 
each graph is obtained from the lower-limit assumption, 
/3=0; the highest curve gives calculated upper-limit values 
based on the use of equations (80) and (81) in calculating 
auL and fiuL] the middle (dashed) curve shows the results 
obtainable by using for a and /3 values determined experi- 
mentally on specimens like those in figure 5. In each case it 
was assumed that /3' = /3, and a' was computed from equa- 
tions (93) and (94). Table I summarizes the upper-limit and 
experimental values of a and /? used for these calculations. 

In general, figures 6 and 7 indicate that the agreement 
between calculation and experiment is within the experi- 
mental scatter, with the calculations based on the values 
Uezp and fiezp giving the best results. 

CONCLUDING REMARKS 

On the basis of an idealization of integrally stiffened plates 
to more uniform plates resembling plywood, formulas have 


TABLE I. 


VALUES OF a, a', /3, AND /S' USED IN THE CALCULATION 
OF THE ELASTIC CONSTANTS FOR COMPARISON WITH 
EXPERIMENTAL MEASUREMENTS OF D*, AND D*, 


hwit-s 

Of exp 

OtUL 



0exp 

5ul 

bwl^s — 9.2 
(a) 

1 

0. 24 

0. 25 

0. 45 

0. 25 

0. 20 

0. 53 

2 

. 17 

. 15 

. 33 

. 24 

. 23 

. 45 

4 

. 12 

. 085 

. 43 

. 31 

. 14 

. 29 

8 

. 004 

. 046 

. 53 

. 43 

. 12 

. 19 

bwibs — 0.4 
(b) 

4 


0. 14 


0. 44 

--- 

0. 14 


» These values (computed from eqs. (74), (77), (80), (81), (93), and 
(94)) were used for calculating constants for all configurations given 
in figures 6 and 7 except those for which 0 — 0° and 0=90° (one-way 
stiffening). 

^ These values (computed from eqs. (74), (77), (81), (87), (93), and 
(94)) were used for calculating constants for configurations of figure 6 
having 0 = 0° and 0 = 90°. 

been derived for the elastic constants of the plates with 
integral ribbing in one or more directions. Two sets of 
elastic-constant formulas have been given, based on two 
different forms of the force-distortion equations. 

The formulas for the elastic constants involve four co- 
efficients a, /9, a', and /3' for each rib which define the effec- 
tiveness of the rib in resisting stretching and bending in its 
transverse direction, horizontal shearing, and twisting. 
Experimental means of determining these coefficients are 
discussed, as are theoretical methods of obtaining values 
corresponding to lower-limit or upper-limit assumptions 
regarding the stiffness of the plate. 

The predictions of the formulas for four of the elastic 
constants are compared with experiment and good correla- 
tion is obtained when experimentally determined values (or, 
in most cases, upper-limit values) of a and /3 are used in the 
formulas for the elastic constants. Despite experimental 
scatter, the calculations and experiments agree, in general, 
both in magnitude and in regard to trends resulting from 
variation in angle of skew of ribbing or in skin thickness. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 26, 1953. 



APPENDIX A 


RELATIONSHIPS BETWEEN NEW AND ORIGINAL ELASTIC CONSTANTS 

The relationships between the new and original elastic constants are as follows: 

Dr 




1 

D. 




1 — 

D 


Di— Z?i(l MxMv) 
Dt=Di{l—nxHv) 


Z>x»=2Z?t 


JS,=- 


Er 


Er 


E, 

^ I— /nf /" 

l-j-jE 


. r ^ r Gu-ura^ t. ^ c^^-y-vCx^ 


Ex—- 


Er 


Er= 


(i^,) <«'->+'A)+o» (t^) (0„+,.cS\ 

E2 

r ^ r Ci2-pirC,, -],r, r C22-JXyCu-\\ 


l+£^; 


Gt, yj GtD^ 

^'‘~l-2Dr,Gryr ^ D^+C^G, 


Cn^Crr ( , A + UzCyr (t 

“ \1 — Mi/*!-/ 

Cil Mx^21 


Cxx 




C2l = IXyCrr f , 

\ j. f*’Xt"y/ N- r *r my 


Dy \ 

fj - 

_ 0\2 

l — l^zfty/ 

^zy 

“Z?i(l — Mi/*») 

Dy \ 

n - 

Gil — HyGll 


K^yX~ 

D 2 ^ 1 M i'} 


(Al) 

(A2) 

(A3) 

(A4) 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 
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C 22 — 




Dy \ 

i — ltzfiy/ 


C22 — l^yCli 
1?2(1 — MiMk) 


(AlO) 


C,= -DryT 



(All) 


=m',+ 

-^)(C„+M + Cyx(Y^^) 

\1 — MiMi// 

■N 

{Cyx’V MxCra;) J+-2Ja: j^(7xx ^ J M M ^ Ml/^w) “1" Gyx ^ ^ ^ ^ ^ (Cyy’\~ M^rCcyJ J 

Ml — 

r C12 — MxC22 

Lz>ill — 

1— 

I+C21 J 

(i-Mxm3 + (i-Mxm3 

C 22 — My^i 2 *l\ 

.Z?2(l— M xM»)J/ 

1+^1-^ On 1 

^ (7ii-MxC'2. -| 

+ C721 [ 

I 1 

^ 3 

fl 

1} J 


(A12) 


M2=/*'»+ 


(Ov+ HyCyy) + m. 6 V)+ O, (j:^) (<?.*+ 
l-Ey (<?*,+ + iOyy+HxCxv)^ 


M 2 — ^ 2 -^ C ^12 

r C\\ — MxCjzi 

LZ>i(l — MxM»)_ 

+C22 

C21 — MyC^U 
_Z?2(1 — MxM»)_ 

} 

1+^2 -^t?12 I 

C\2 — MxGj 2 "1 

_Z?l(l — /lx/»»)J 

+C22 

C 22 iXyC\2 

_Z>2(1 — 

r> J 


(A13) 


APPENDIX B 

DERIVATION OF FORMULAS FOR ELASTIC CONSTANTS 


The basic assumptions of the analysis have already been 
described. In the derivations that follow, where the word 
‘‘rib^’ is used, it means one of the substitute sheets, depend- 
ing on which property of the rib is under consideration. 
Separate derivations are given for the constants associated 
with bending and stretching and those associated with twist- 
ing and shear. 


CONSTANTS ASSOCIATED WITH BENDING AND STRETCHING 

In the derivation of the formulas for the elastic constants 
•associated with bending and stretching, an element of the 
integrally stiffened plate will be considered; the element has 

the average prescribed curvatures and and the strains 

€x (measured in some arbitrary plane which will be referred 
to as plane I) and €y (measured in some other arbitrary plane 
which will be referred to as plane II). The development of 
these prescribed deformations requires the application of 
moments of intensity Mx and My and forces of intensity Nx 
(acting in plane I) and Ny (acting in plane II). These mo- 
ments and forces and the locations of planes I and II are 
shown in figure 8 . 

If the strains are assumed to vary linearly through the 
thickness of the element, two horizontal planes can be found 

^in terms of ej, and in which the x-wise strain 

and y-wise strain, respectively, are zero. These planes are 
indicated in figure 9. 

Strains of components of plate. — The longitudinal exten- 
sional strains of the ribs measured at their cross-sectional 
centroids can be written in terms of the curvatures and the 
distance between the rib centroids and the planes of zero 
extensional strains. The strains of the x-wise, y-wise, and 


X 





y 



Figure 8. — Forces and moments considered for analysis of bending 

and stretching. 


skewed ribs are, respectively, 



(Bl) 

€w — ^ 0 

(B 2 ) 


(B3) 


where the subscript L denotes longitudinal direction of a rib, 
the subscript x the a:-wise rib, the subscript y the i/-wise rib, 
and the subscript s the skew rib. The distances A 3 , Ara, Ai, 
and ki are shown in figure 9. 

The transverse strains of the ribs are as follows: 


f, rrv 

= — ikn — - — S- 

(B4) 

^W^=-{h-C^ytI) ^ 

(B5) 


— sin^ e—{k2—ajl) cos^ B (B6) 


The extensional strains of the sheet jnidplane in terms of 
the curvatures are 




(B7) 


^S,——k2 ^ 




(B 8 ) 


The curvatures and of the element are also the cur- 



' 1 

z^w 


y 



• Rib centroids 

Figure 9. — Dimensions for analysis of bending and stretching. 
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vatures of the ai-wise and y-wise ribs, respectively. The 
curvature of the skew ribs is 


^=_cos*0+^sm»e 


(B9) 


and At 2 from equations (B13) and (14) gives 


(B19) 


The horizontal shear strain in one of the skew ribs, relative 
to the longitudinal and transverse directions of the rib, can 
be written in terms of the ar-wise and y-wise strains at the 
same level, which in turn are determined by the a:-wise and 
y-wise curvatures; thus, 

yw=2 1^— a',H) ^+(At 2 — a',H) sin 6 cos 6 (BlO) 

The x-wise and y-wise ribs have no shear strain. 

Expressions for the dimensions hz, ku At 2 » and kz * — In 
the derivation of equations (Bl) to (B8) and of equation 
(BIO), the assumption was made that the strains varied 
linearly from the planes of zero strain. On the basis of the 
same assumption, expressions are written for the strains in 
planes I and II — the planes in which and Nj, act and in 
which €x and ty are measured. These expressions are 

(Bll) 
(B12) 


(B13) 




e^= — {ki—knH) 


from which 


bf 


fh=kiH- 


b^w 

bl? 


hz={k^-h)H^^- (B20) 


(B21) 

W 


ki— (kwy—kn) H+ (B22) 

bf 

Evaluation of strain energy. — The total strain energy of 
the element of the integrally stiffened plate can be written as 
the sum of the strain energies of its component parts; thus 

-J^ ^w,^^EAw,ds+-^^^ ^^zrE^zAwjix-\- 

1 1 r&*8«oe 

2Jo ^^»TE^>^w,ds-\- 





h=kuH-^ 

bf 


(B14) 


By geometry the dimensions hi, h^, ki, and ks may be 
written 


h\ — ky^Il — 


(B15) 




hz — kw^H — ^2 (B16) 

k\=kj^H — A*2 (B17) 

kz—kwyH kz (B 18 ) 

where kiryH, kwH locate the centroidal axes of the 

ribs from the center line of the sheet. Substituting for 


In equation (B23) the first three terms give the energy of 
extension of the ribs in their longitudinal directions, the 
second three terms the energy of extension of the ribs in their 
transverse directions, the seventh term the energy associated 
with the shearing of the ribs, and the eighth term the energy 
of extension of the skin. The next three terms give the 
energy of bending of the ribs, and the final term gives the 
energy of bending of the skin. 

Carrying out the integrations of equation (B23), dividing 
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by hzh^ to reduce the result to strain energy per unit area, and substituting the previously derived expressions for the 
distortions ew,, and so forth gives 

bzb, 




^i) H cos^ ^-4" 




/3jsin^0cos^5— /3',Y:^sin“ ^ cos^ ®)J ArO H+^y"^iay—k^ H+^~^(Jcw—ki)H cos* 6+ 
/3, (a,—ki) H sin‘ d+fi', {a',—ki) H ® ** (.hiH)+^^ (kw,—ku) H sin^ B cos* 0-f 

0, {a>—kii) H sin* e cos* B—fi'z (a',—kn) h(^ sin* 6 cos* <s+/3x ~^+~p^+ 

-^^sin^ 0+/3, cos^ $+i8'jY:^®*’^^ ^ t„*-h2-|^ (kjH)-\ — ^ (kw,—kj) H sin^ 6 cos^ 6+ fi, (as—kt) H sin^ 6 cos* 6— 

/3'. ia'-k^) h(y|yu ^ ^)]} S^+^T^* H+ 

Atii) H sin^ e+^, (a,—ku) H cos‘ B+^\ {a' ,—kn) H ^ 1. 12 (1— 

-Zii^ iw 1 -^4. w — w .i4. IV r — 

X"'^x ®+I^7* (^H^^-A:,)*^/*-f /3, (ay-kr)m^+-^ [^(A:^ -Ari)*i?* cos‘ B+fi, {a-hfH^ sin^ B+ 

/3'. (a%— A ti)*//* (^y:|^ ® ~*~^^12 (1— 7*) ^ (A:iA:ni/^+ 

-4,r r - 

-y^ I A;i) (^tr,— A tii) sin* B cos* 04-/3, (a,— Aij (a,— Arn) sin* 0 cos* 0— 

r, (a'.-Arx) (a'.-Ar,:) /^* (y|^ sin* 0 cos* 0)]}^ ^4- {127^ ^^'+X“+X 
-^d-w — x4ur r — 

/3x («X-A:n)*fl'*4-JY-" ik^-knfH^+-^ I (Ar,r -Arnl^i?' sin^04-/3. ioL-ki^fW cos* 0+ 

7. (a',-A:xi)*H* (y-|^ sin* 0 cos* «)]}(0j) (B24) 

where the identities ^ ^ have been substituted to simplify the expressions. 

Invoking the principle of virtual displacements by differentiating the energy expression (B24) with respect to each 


j^w r — 
T-* 
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of the strains and curvatures and dividing by EH or EW gives the following expressions for the forces and moments: 


W 1 

dex EH~EH 


(cos‘ e+p. sin‘ 0+^’, sin^ $ cos^ o)] **+[1372 ^ 

fi, sin* 0 cos* 0—^', ^ ^ H — W~~ («»“^i)+ 

cos‘ e+|8, (a.—kj) sin* 0+^', (a',-kj) sin* 0 cos* ^ *h+ 

^ii) sin* 0 cos* fl+/3, {a,—kii) sin* 0 cos* 0—fi', (a',—kn) sin* 0 cos* ^ 


d*w 

w 


(B25) 


i>V 1 Ny 

d«„ EH~EH 


(^sin* 0 cos* 0+p, sin* 0 cos* 0—p', sin* 0 cos* " ~l — 

^sin^ 9+/3, cos‘ 0+fi', sin* 0 cos* *^+ ^ ki-\ — Ati) sin* 0 cos* 0+ 


P, {a,—ki) sin* 0 cos* 0—P', {a',—kj) sin* 0 cos* H («i— Arn)+ 

^ Ocw.—kn)+^^~ j^(*n/.— Aril) sin* fl+/3, (a,— Ar„) cos* 0+P', (a',—ku) sin* 0 cos* ^ (®26) 


Aw^ 

~Tl 


dF' 1 Mx 

^d^wEH^~ EW 
^8? 


1 h ,. , ‘^wjh, ,r 


^jF/ H 


( kw —^ i)-\-^y ^ { oiy — ki )-\ Ati) cos* 6-\-Pg ( ag — ki ) sin* 6 + 


P', {a',-ki) sin* 0 cos* «x+ ^ j^(^„.,-A:i) sin* 0 cos* 0+p, (a.-kj) sin* 0 cos* 0— 

P', sin* 0 cos* 0^^ (y+ -[i2 (1-M^ (h) '^+1^* (*ir,-A:i)*+ 

A:i)* cos* 0+p, (a,— A ti)* sin* 0+P', (a',— A tj)* ^ 

0 cos* 0 ^ kikn+ ^-^- A:,) (kw-ku) sin* ^ cos* 0+ 




/ M /t«Y , 


p, (ot,—ki) (a,—kii) sin* 0 cos* 0—p', (a',~ki) {a',—kn) ® 


d*w 

dy* 


(B27) 
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bV' J. Mr 

EtP 


= ^ ^ 11 +-^- ^Ii) sin* e cos* e+fi,{a,—kii) sin* $ cos* e—p',(a',—kn) sin* $ cos* «*+ 

{l^* h sin^ 0+/3. («.-*„) cos^ o-h 

/3'. («'.-Ar,0 sin* e cos* { j 2 (i^ (sj ^ ''+1^ I 

^ j^(itir,— i^i) {kw,—kii) sin* 6 cos* d+ft(a,— lri)(a,— ^ii) sin* $ cos* ^— 


/3'.(a'» *i) («'. *n) sin* 0 cos* H 12(1-^*) (ff) 


1 .:.4 fl . _J_ £ 


^ ) +-W+-W *^+^* ^ 


1-u* H 


^ Awjby 7 \2 I 

Pz JJ (^z ^II) ”1" u 


/r 

H 


{kw^—kiiy+'^^^ ^(*ir,— *ii)* sin< ^+ 


^.{ct.—kuy cos‘ 0+fi',(c 




The equations for N^y Nj,y and M„ (eqs. (B25) to 
(B28)) can be written as 

^=A,€,+A.^+A,(k,-h)H ^+A.(l-k^,)^ 0 (B29) 
■^=A.e,+A,^+A.(k-h)H^+Ay(ky-k^^)H^ (B30) 

[l.+A.gi.-h){k.-ki^)]H^ (B31) 

^II) ^Z~\~ -Ay(ity ^II) ^V~\~ 

[T.+A.(ic.-h)(k-ku)]H g+ 

(B32) 

where Axy A, and so forth are given in equations (43) to (54). 

In order to identify the desired elastic constants associ- 
ated wdth extension and bending, the foregoing force- 
distortion relationships, equations (B29) to (B32), need only 
to be put into the form of equations (1), (2), (4), and (5) or 
equations (7), (8), (10), and (11). 

CONSTANTS ASSOCIATED WITH TWISTING AND SHEARING 

The derivation of the formulas for the elastic constants 
associated with twisting and shearing is a parallel one to that 


for the bendmg and stretching constants. 

An element of the integrally stiffened plate which has the 

average prescribed twist shear strain 7xy (measured 

in some arbitrary plane which is referred to as plane III) is 
considered. These prescribed deformations can be effected 
by the application of twisting moments of intensity M^y and 
shearing forces of intensity Nxy (acting in plane III) to the 
element. (See fig. 10.) 

If the horizontal shear strain is assumed to vary linearly 
through the thickness, the horizontal plane can be found 



2 Plane m 


Figure 10. — Shears and moments considered for analysis of twisting 

and shearing. 


22 


REPORT 1195 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


0 


in terms of 




whi 


and ) which has zero shear strain. 


da: dy 

This plane is shown in figure 1 1 . 

Strains of components of plate. — The extensional strains 
of the longitudinal, transverse, and one of the skew ribs in 
their longitudinal directions at their centroids are 

(B33) 

(B34) 

tw =±h\ sin 2$ (B35) 

•L i>xdy 

The transverse strains of the ribs are 

=0 (B36) 

Xj, 

(B37) 

± {h' 2 —oisH) ^ sin 26 (B38) 

The extensional strains of the sheet are 

€5,=0 (B39) 

65^=0 (B40) 

The twist causes bending of the diagonal ribs. The 

curvature of one of these ribs is given by 

The curvatures of the longitudinal and transverse ribs are 


zero. The shear strain in the skin middle surface is given by 


7s=— 2A'2 


dx by 


(B42) 


The magnitude of the shear strain of the diajgonal ribs is 
given by 

Tpr,=2(A'2— COS 26 


(B43) 


The shear strain of the a:-wise and y-wise ribs is given by 

b^w 


— — 2{hf 2 cx! xH) 
yw =— 2(A'2— 


bx by 

b^w 
bx by 


(B44) 


(B45) 


Expressions for the dimensions h\ and — The following 

expressions can be written for the strains in plane III, in 
which Nxy acts and in which ^xy is measured (see fig. 10): 


yxv — — 2{N2 — kiiiH) 


b'^w 


from which 


bx by 


hf 2 — kiiiH — - 


2 b’^w 


By geometry 


bx by 

h\=k,vH-h'2 

Substituting for h '2 from equation (B47) gives 


(B46) 

(B47) 

(B48) 

(B49) 


di/ 



z, tv 



Evaluation of strain energy. — The total strain energy can 
be written as 


A. w.ds-\~ 


. 1 rs.secS j , 1 r*-*- 

“2 Jo Jo 

2 J^ '^w^G^'xAwdx-\--^ J^ 




i r- r 

^ Jo Jo 


y sGtf^dxdy-\- 




In equation (B50) the first term gives the energy of 
extension of the skewed ribs in their longitudinal directions; 
the second term, the energy of extension of the skewed ribs 
in their transverse directions; the next three terms, the 
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energy of shearing of the ribs; and the sixth term, the energy of shearing of the skin. The next term represents the energy of 
twisting of the skin, and the last term gives the energy of bending of the skew ribs. 

Carrying out the integrations of equation (B50), dividing by hxhy^ substituting previously derived expressions, and so 
forth, gives 


U 


hxh\ 








^ /3'x^ (oc'-knt)H+^^^ (liw-h^^)H sin^ d cos^ B+ 

fi.(a,-hu)H 8in» 6 cos^ e+^'.{a'.-km)H [ 2 ^ cos^ 2«]}) ^+G(1^ 


4 ^ sin* e cos* ^'x («'x-^„i)^H ^+ &'v ^ 

4 ^ { (kv,-k^^,YW sin* e cos* e+p,{a,-kj„YW sin* 6 cos* 0+/3',(a',-Ar„,)*H* [^!p^ cos* 2d] ) (^1^)'] 

Differentiating the energy expression (B51) with respect to each of the distortions and dividing by EH or £H* gives the 
following expressions for the forces and moments: 

HU' 1 ^ 

HyiryEH EH 

^ i, 2(1+^ ■^2(1+7) ~HT'^~il~ ^ ^ 

2(1+7 ]“ Txv+ 2 (2(1+7 ^*«i+2(r+^ 

— 'Yj ' ■( (kw,—kiii)sin^ d cos* 6+^, (a,— km) sin* d cos* d+^',(aa— Arm) [^ 2(i+ju) ^ Hx'dy (B52) 


HU' 1 _2 


d*w EH^~ EH^ 


bx by 


(a'»— Arm)-! — A:„i) sin* d cos* d+ 


_o{ —1 7. I 1 ot ( > 1. 1 f,, ^w+ix , 

\2(1+M)fl’*'"'''2(l+M)^* H *'"^+2(l+Mr'' H 

/3,(a.— A tiij) sin* d cos* d+^'.(a's— A-m) |^ 2 ~ (i +7 } ^xy+( ^j^^j (^^ +4 sin*d cos* d- 

2 2 -^4. nr / hx .id. nr /&• ^ — 

14 ^ S ^ — ^iii)^ sin^ 0 cos^ 6-\r 

/3,(a,— A:i„)* sin* d cos* d+/3',(a',— Arm)* [ ^2([+7 bxby 


^ H(A;„,)*+ 


1+mA?' 


(B53) 
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The equations for and M*„ (eqs. (B52) and (B53)) can be written as 

^=A^y^^2Ajh-hu)H ^ (B54) 

2 ’^^—^Azfikxy kiii)‘yxv~\‘\Ixy~\'^Axy{kxv (B55) 

where and hy are given in equations (46), (50), and (54), respectively. 

Equations (B54) and (B55) may readily be put into the form of equations (6) and (3) or (12) and (9) to yield either 
the original or the new elastic constants, respectively. 
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